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1. Introduction 

Study of the p meson decay is a significant step for understanding the dynamical aspects of 
hadron interactions with lattice QCD. In the early stage of studies toward this direction the transi- 
tion amplitude (p\nn) extracted from the time behavior of the correlation function (n(t)n(t)p(0)} 
has been used to estimate the decay width, assuming that the hadron interaction is small ^ ||,|3|, |J. 

A more realistic approach is an estimation of the decay width from the P-wave scattering phase 
shift for the isospin 1=1 two-pion system. The finite size formula presented by Liischer in the 
center of mass frame |Q] or its extension to non-zero total momentum frame by Rummukainen and 
Gottlieb [J7|] is employed for an estimation of the phase shift. The first study of this approach was 
carried by CP-PACS Collaboration using Nf = 2 full QCD configurations {m n = 330 MeV, a = 
0.21 fm, La = 2.5 fm) [0]. After this work, two studies were reported with Nf = 2 configurations 
on finer lattices, one by QCD-SF Collaboration (m„ = 240 - 810 MeV , a = 0.072 - 0.084 fm) [§] 
and the other by ETMC Collaboration (w% = 391 MeV, a = 0.086 fm, La = 2.1 fm) fl9|. 

In the present work we extend these studies by employing Nf = 2 + 1 full QCD configurations 
and working on a larger lattice volume. Our calculations are carried out on a subset of configura- 
tions previously generated by PACS-CS Collaboration with the Iwasaki gauge action and nonper- 
turbatively 0(a) -improved Wilson fermion action at j3 = 1.9. on a 32 x 64 lattice [10]. The subset 



corresponds to the hopping parameters K u d = 0. 13754 for the up and down quark and K s = 0. 13640 
for the strange quark. The parameters determined from the spectrum analysis for this subset are 
m % = 410 MeV {m n /m p = 0.46), a = 0.091 fm and La = 2.91 fm. All calculations are carried out 
on the PACS-CS computer at Center for Computational Sciences, University of Tsukuba. We note 
that ETMC [11] and BMW JI^ ] Collaborations reported their preliminary results at Lattice 2010. 



2. Method 

We consider the center of mass frame (CM) and the non-zero total momentum frame (the 
moving frame (MF)) with the total momentum p = (2n/L)ej. In these frames the ground (n = 1) 
and the first exited states (n = 2) with spin 7 = 1 and isospin 1=1, ignoring hadron interactions, 
are given by 

n = 1 n = 2 

OM n(ej)n{-Bj) [1.3] 

Pi 7(0,0,1) n(\,0, 1)^-1, 0,0) , 7r(0,l,l)7r(0,-l,0) [1.4] ' 

galQAl) 7r(0,0,l)7r(0,0,0) [1.02] 

(2.1) 

where p is the total momentum, g is the rotational group on the lattice and T is the irreducible 
representation of the group. The vectors in parentheses after % and p refer to the momenta of the 
pion and the p meson in units of 2%jL. The numbers in square brackets of the two-pion states are 
values of ^/s/m p on our full QCD configurations. In the present work we calculate the scattering 
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phase shifts of the states marked by under-bar in (2.1). The finite size formulas for these states are 
given in Refs. [g|]. 

For the and the E~ representations the energy of the ground state is much smaller than 
that of the exited state as one can see from the value of y/s/m p in ( |2. 1[ ). Thus the energy of these 
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states can be extracted by a single exponential fit for the time correlation functions of the p meson. 
We use the local p meson operator for the sink and a smearing operator for the source as discussed 
later. 

For the Aj representation the energy of the first exited state is close to that of the ground state. 
Thus we use the variational method [IT3|] with a matrix of the time correlation function, 



G(t) 



(2.2) 



(0^)0(4)10) (01^(0^(010) J ' 

We extract the energy by a single exponential fit for the two eigenvalues X\ (t) and Xz(t) of the ma- 
trix M(t) = G{t)G~\t R ) with some reference time assuming that the lower two states dominate 



the correlation function. In ( |22| ), p^(t) is the local operator for the neutral p meson at time t with 
momentum p and the polarization vector parallel to the momentum. Cl(t) is an operator for the two 
pions with the momentum and p = (27i/L)e3, 



(p,t)-n (0,ri)w + (p,f))xe 



jn n -{t] -t) 



(2.3) 



where n(p,t) is the local pion operator with momentum p at time t. The times slice of the pion 
with zero momentum is fixed at t\ » t, and an exponential time factor is introduced, so that 
the operator has the same time behavior as that of the usual Heisenberg operator for t\ » t, ie., 
(0|£2 + (f) = (0|ft t (0)exp(-H-f)- 

Two operators pj(4) and £l(t s ) are used for the source in ( |2.2| ), which are given by 

(2.4) 
(2.5) 



Q(f,) --^(;; I) / (p./,) - k (0./, i7r + (p./ v ) 



■D(z,t s )ftD(z,t s ))e 



zer V 2 v 

The operator U (z,t s ) (D(z,t s )) is a smearing operator for the up (down) quark given by U (z,t s ) = 
^ x u(z + x, t s ) ■ F(\x — z\ ), where u(x, t s ) is the quark operator at position x and time t s . We use the 
same smearing function F(x) as in Ref. [[K]]. This operator is used after fixing gauge configurations 
to the Coulomb gauge. In ( |2~3| ) a summation over z is taken to reduce contaminations from the 
states with different total momenta and T = { z | z = (L/2) • (n\,ri2,m) , rij = or 1 } is chosen in 
the present work. The smearing operator ( |L5| ) is also used to extract the energy of the ground state 
for the Tj~ and the E~ representations setting the momentum to p = and (2n /Vjej. 



The quark contractions of G(t) in ( |2.2[ ) are given by 
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(2.6) 
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where the time runs upward in the diagrams. The vertices refer to the pion or the p meson operator 
with the momentum at the time specified in the diagrams. The p meson operators at t s are the 
smearing operators, and the other is the local operator. 

In order to calculate the quark contractions in {!£), we use the source method and the stochas- 
tic noise method as in the previous work by CP-PACS Collaboration [0]. We introduce a U{\) 
noise £,j(x) which satisfies 



Nr 



£y(x)^(y) = 5 3 (x-y) for N R 

7=1 



(2.7) 



where Nr is the number of noises taken to be 10 in the present work. We calculate the following 
four types of quark propagators : 



Q AB (x,t\q,t s ,^j) =£(D-V(x,?;y^) • e*^(y) 
y 

W AB {x,t\k,t a \q,t s ,$j) =££(D- 1 ) AC (x,f;z,f a )- [e*^ Q(z,t a \q,t s ,$j) 
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CB 



Q AB (x,t\z,t s ) = 1 £{D l ) AB (x,t;y,t s )- 



niy-z| 



W AB {x,t\kf a \z,t s ) =££(£> l ) AC {xf,yf a ) 
y c 



e ik - y 75 Q(y,t a \z,t s ) 



CB 



(2.8) 
(2.9) 
(2.10) 
(2.11) 



where A and B refer to color and spin indices, and F(x) is the smearing function. The square 
bracket is used as the source for the inversion of the Dirac operator D. 

The first term of the kk — > kk component of (03) can be calculated by introducing an another 



t/(l) noise T] ; (x) having the same property as ^y(x) in < p?7) ), 

EEe-WfitfctilO,^) e(x,? 1 |0, fs ,^))-e'"-("-').(G t (y,?|p,?. v ,T 7; -) QMvM)) 

j *,y 

(2.12) 

where the bracket means trace for the color and the spin indices. The exponential time factor 
comes from the definition of the operator of the two pions in (2.3). The second term is given by 
exchanging the momentum and the time slice of the sink in (2.12). The 3rd to 6th terms can be 
calculated by 

g££U(0 =££e- ! > x .(wHx,?|0, ?1 |0,^,^)W(x, ? |0,? ? |p,^,^))-e m -^- ( ), 



4T-U(0 = £!>-*■*.( wHx,t\0, tl \0,t s ,Sj) W(x,t\p,t,\Q,t„tj) )-e 



G [ nnLn(t) = £JV** ' ( W(x,t\0,h \0,t s , $j) W\x,t%t s \ - \>,t s , fy) ) ■ e 



J * 



jn n -(t]_-t) 



= lie-*'" • ( W(x,t\0,h\0,t s ^j) W\x,t\-p,t s \Q,t s ^j) Y e m ^^ .(2.13) 
j x 



The two terms of KK — > p of ( |2.6| ) can be similarly obtained by 

>St \ P [t) = 

j 



G l £l p (t) = £Ee-*-*-(w t (x,t|-p,t,|0,r„^) (7573) <2M<U,^ 
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gKp(0 = ££e-* x - ( ^MO,^;) (7573) W(V|p,**|<M„&) ) . (2.14) 



,[2nd] 

•pv; — 

We can calculate the two terms of p — >• kk of ( |2.6[ ) by 

4^U0 = £e*"£e-*-«.(S t (x,^,f,)^(x,f|0 > f 1 | Z) O (737s) )-e m - (;i - () , 
zer x 

Gj^UO = £e i »"£e-*- x .(w + (x,/|0,fiKOfi(*.^.^) (737s) )-e m "-("-') . (2.15) 
zer x N ' 

The p — >• p component is given by the <2-type propagators as the usual two point function. 

The quark propagators are calculated with the Dirichlet boundary condition imposed in the 
time direction and the source operators are set at t s = 12 to avoid effects from the temporal bound- 



ary. We put the zero momentum pion introduced in (2.3) at t\ = 42. We calculate the £2-type propa 



gators (2.8) for four combinations of q and U(l) noise : (q, noise) = {(0, £), (0, T}), (p, E,), (— p, £)}. 
The W-type propagators ( |2.9[ ) are calculated for five combinations of k, t a and q : (k,f a |q) = 
{(p,?s|0), (— p,ft|0), (0,ts\p), {0,ts\ — p),(0,?i|0)}, using the same U(\) noise t, in common. We 



calculate the 2-type propagator ( |2.10| ) and the VK-type propagator ( 2.1l[ ) with (k,f a ) = (0, ?i ) for 



the set z E r. Thus we calculate (4 + 5)xl0 + (l + l)x8=106 quark propagators for each config- 
uration. The total number of the configurations analyzed every 10 trajectories is 440. We estimate 
the statistical error by the jackknife method with bins of 400 trajectories. 



3. Results 

In the left panel of Fig. [I] we show the real part of the diagonal components (nn — > nn and 
p — > p) and the imaginary part of the off-diagonal components (nn — > p, p — > Tin) of G(t). The 
other real or imaginary part of the components vanish from P and CP symmetry. We calculate 
the two eigenvalues X n (t) (n = 1,2) for the matrix M(t) = G(t)G~ l (tn) with the reference time 
tR = 23. In the right panel of Fig. [j] we plot the eigenvalues normalized by the correlation function 
of the two free pionsN(f) = (0|^(-p,f)^(p,^)|0)(0|^(0,f)^(0,4)|0). Thus the slope of the figure 
corresponds to the energy difference with respect to the energy of the two free pions. We observe 
that the energy difference is negative for X\ (t) and positive for ^(t). This means that the two-pion 
scatting phase shift is positive for the ground state and negative for the first exited state of the A2 1 
representation. 

We extract the energies for both states by a single exponential fitting of the eigenvalues X n (t) 
(n = 1,2) for the time range t = 24 — 33. For the T~[ and the E~ representations the energy of the 
ground state is obtained from the time correlation functions of the p meson as explained before. 
Converting the energies on each frame to the invariant masses and substituting them into the 
finite size formulas, we obtain the scattering phase shifts plotted in Fig. ||. In the figure we show 
(k 3 / tan8(k))/y/s as function of (y/s) 2 in unit of the lattice cutoff, where k = y/s/A — m^ is the 
scattering momentum. 

In order to estimate the p meson decay width at the physical quark mass we parametrize the 
scattering phase shift with the effective p — )■ nn coupling constant g p7[K by 

k 3 , r- 67T . 9 

—— yf S =—. {M\-S , (3.1) 

tan 5 (A:) «L 
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Figure 1: Left panel : Real part of the diagonal components (nn — > nn and p — >• p) and imaginary part of 
the off-diagonal components (jttt — > p and p — > jttt) of the time correlation function Git) ; Right panel : 
Eigenvalues X\ (t) and Xi{t) normalized by the correlation function of the two free pions N(t). 
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Figure 2: (A: 3 / tanS(k))/y/s as function of (\/s) 2 in unit of the lattice cutoff. Dotted lines refer to the finite 
size formulas for each representation. 



where Mr is the resonance mass. The coupling g pK7[ generally depends on the quark mass, but our 
present data at a single quark mass do not provide this information. Here we assume that the mass 



dependence is weak and try to estimate g pK7[ and Mr by fitting our results with (p.l[). We estimate 
the p meson decay width at the physical quark mass from the formula F p = g 2 p nn x 4.128 MeV. 
Our results of the fitting are given by 

oiMr = 0.4064 ± 0.046 , g p%% = 5.24±0.51 , r p = 1 13 ± 22 MeV . (3.2) 

In Fig. H we indicate the result for the resonance mass Mr and draw the fitting curve by a solid line. 
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Our result of the p meson decay width T p at the physical quark mass is smaller than the 
experimental value (150 MeV). Possible reasons of the discrepancy are the mass dependence of 
the effective coupling constant g pnn and the finite lattice spacing effects. We will study these issues 
in next investigations. 

4. Summary 

In the present work we estimated the p meson decay width from the Nf = 2 + 1 full QCD 
configurations generated by PACS-CS Collaboration. The decay width is estimated from the P- 
wave scattering phase shift for 7 = 1 two-pion system. We used the effective p — > nn coupling 
constant g pn7[ to extrapolate from our simulation point m n = 410 MeV to the physical point m n = 
135 MeV, assuming that g pK7l does not depend on the quark mass. The decay width may be 
estimated directly from the energy dependence of the phase shift data assuming the Breit-Wigner 
resonance formula, if the simulations are made close to the physical quark mass and we have data 
at several values of energy near the resonance mass. We leave these issues to studies in the future. 

This work is supported in part by Grants-in-Aid of the Ministry of Education (Nos. 
20340047, 20105001, 20105003 , 20740139 , 20540248 , 21340049 , 22244018, 20105002 , 
22105501, 22740138 , 10143538 , 21105501 20105005 ). The numerical calculations have been 
carried out on PACS-CS at Center for Computational Sciences, University of Tsukuba. 
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